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Abstract 

First theoretical analysis of the nuclear spin conversion in ethylene molecules 
(^■^00114) has been performed. The conversion rate was found equal ~ 3 • 
10~^ s~^/Torr, which is in qualitative agreement with the recently obtained 
experimental value. It was shown that the ortho-para mixing in ^'^CCH4 is 
dominated by the spin-rotation coupling. Mixing of only two pairs of ortho- 
para levels were found to contribute significantly to the spin conversion. 
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I. INTRODUCTION 



Nuclear spin isomers of molecules were discovered nearly 70 years ago when ortho and 
para hydrogen isomers were separated for the first time. Although it was realized already at 
the time of this discovery that many other symmetrical molecules have nuclear spin isomers 
too, until recently almost nothing was known about isomers of molecules heavier than hy- 
drogen. This happened because of the lack of practical separation methods. The separation 
method based on deep cooling was applicable to hydrogen (and deuterium) but failed in 
the case of heavier molecules. Presently a few new methods for spin isomer separation have 
been proposed and successfully tested which has advanced this field significantly. Yet we 
are at the very early stage of this research. The list of molecules for which the nuclear spin 
isomers have been separated is rather short: H2 [0], CH3F §, H2O @], CH2O §, Lis §, 
i, and C2H4 i. 

The CH3F molecules occupies a special place in this list because it is the only poly- 
atomic molecule for which the isomer conversion mechanism has been identified. The CH3F 
spin isomers appeared to be extremely stable surviving ~ 10^ gas collisions. Nevertheless, 
the isomer conversion was found to be based on ortho-para state mixing induced by tiny 
intramolecular hyperfine interactions and interruption of this mixing by collisions. This 
specific type of relaxation was proposed to refer to as quantum relaxation. It is essential for 
the physics of molecular spin isomers to understand which processes are responsible for the 
spin conversion in other molecules. 

Recently first separation of nuclear spin isomers of ethylene molecules (^^CCH4) has 
been performed 0. The conversion rate has been determined and it was shown that the 
rate increases proportional to the gas pressure, 7/P ~ 5 • 10^"^ s^^/Torr. An understanding 
of mechanism of the conversion in ethylene needs theoretical investigation of the process. 

In the present paper we perform first theoretical analysis of the spin conversion in ethy- 
lene molecules. The purpose is to verify the consistency of the experimental data on isomer 
conversion in ethylene with the conversion by quantum relaxation. We consider the same 
isotope species, ^^CCH4, for which the experiment was done 0. An essential difference of 
the present study from the well-understood case of CII3F molecules is that now we have to 
work with an asymmetric-top molecule. This appeared to be much more complicated than 
the case of symmetric-top molecules. We would like to point out that spin conversion in 
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asymmetric tops were considered theoretically previously for the CH2O and H2O molecules 
[|]- 



II. ISOMER CONVERSION BY QUANTUM RELAXATION 

Although the spin conversion in molecules by quantum relaxation is explained in other 
publications we will describe briefly its essence here for the convenience of the reader. At 
room temperature almost all molecules are situated in their ground electronic and vibrational 
states. Suppose that these states are separated into two subspaces which are the nuclear spin 
ortho and para states, as is shown in Fig. |l] for the particular case of the ^^CCH4 molecules. 
Note that some molecules may have more than two nuclear spin isomer forms, e.g., methane, 
or normal ethylene. On the other hand, number of molecules, e.g., H2, CH3F, and ^^CCH4 
too, have just two type of spin species. 

The relaxation process which we are going to consider has two main ingredients. First, 
the ortho and para quantum states of the test molecules are not completely independent. 
There is small intramolecular perturbation, V, which is able to mix the ortho and para 
states. This perturbation can mix, in general, many ortho-para level pairs. But in Fig. ^ 
the mixing is present just for one pair of states as will be the case for ^^CCH4, see below. 
Second, the test molecules are embedded into an environment which is able to induce fast 
relaxation inside the ortho, or para subspace, but is not able to produce direct transitions 
between these subspaces. This implies that the relevant cross-section, a{ortho\para) = 0. 
The relaxation by direct transitions is the process opposite to the quantum relaxation. 

The isomer conversion by quantum relaxation consists in the following. Suppose that at 
the instant t = a. test molecule is placed into the ortho subspace. Due to collisions with 
surrounding particles, the test molecule starts to perform fast migration along rotational 
states inside the ortho subspace. This is the familiar rotational relaxation. This running 
up and down along the ladder of the ortho states continues until the molecule jumps in the 
state a which is mixed by an intramolecular perturbation with the energetically close para 
state a'. During the free flight after this collision, para state a' will be admixed to the ortho 
state a. Consequently, the next collision can transfer the molecule in other para states and 
thus localizes it inside the para subspace. Such a mechanism of spin isomers conversion was 
proposed in the theoretical work by Curl et al (see also [^). 
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The quantum relaxation of spin isomers can be quantitatively described in the framework 
of the kinetic equation for density matrix . One needs to split the molecular Hamiltonian 
into two parts 

H = Ho + V, (1) 

where the main part of the Hamiltonian, Hq, has pure ortho and para states as the eigen- 
states; the perturbation V mixes the ortho and para states. If at initial instant the nonequi- 
librium concentration of, say, ortho molecules Spo{t = 0) was created, the system will relax 
then exponentially, Spo(t) = 6po{0)e~'^^, with the rate 

7= E ^^^iW,{a') + Wo{a)); F{a'\a) ^ K'a\'. (2) 

The sets of quantum numbers a' = {a', u'} and a = {a, u} consist of the degenerate quantum 
numbers u', v and the quantum numbers a', a which determine the energy of the states. In 
(ID r is the decay rate of the off-diagonal density matrix element pa'a («' £ para; a G ortho) 
assumed here to be equal for all ortho-para level pairs; Ua'a is the gap between the states 
a' and a; Wp{a') and Wo{(y) are the Boltzmann factors of the corresponding states. For 
the following it is convenient to introduce the strength of mixing, F(a'|a), which sums the 
intramolecular couplings over all degenerate states. 



III. ROTATIONAL STATES OF ETHYLENE 



Ethylene is a popular object in high resolution infrared spectroscopy. This fortunate 
circumstance made available rather accurate data on molecular parameters and position of 
molecular rotational levels. The molecular structure, numbering of hydrogen atoms and 
orientation of molecular system of coordinates are presented in Fig. |^. ^^CCIl4 is a plane, 
prolate, nearly symmetric top having the symmetry point group The characters of the 
group operations and the irreducible representations are given in the Table 1. We give in the 
Table 1 also two isomorphic groups, molecular symmetry group C2v{M) and the point 
group D2. The bond lengths and angles which were used in our calculation are as follows: 
rcH = 1-087 A, rcc = 1-339 A, the angle aiyciy=117.4° JIT]. 



The fact that the ethylene molecule is an asymmetric top complicates in two aspect the 
theoretical analysis of the isomer conversion. First, the energy levels and wave functions 



for asymmetric tops can be found only numerically but not analytically as is possible for 
symmetric-top molecules. Second, the rotational quantum number k (projection of molecular 
angular momentum on the molecular symmetry axis) becomes an approximate quantum 
number in asymmetric tops. Consequently, one should calculate numerically more ortho- 
para matrix elements than in the conversion of symmetric-top molecules, see, e.g., M. 



Rotational states of CCH4 can be determined using the Hamiltonian of Watson 



and the set of molecular parameters from Refs. p!3|JT^ . The rotational Hamiltonian up to 
sextic order terms has the form 



1(5 + C)J' + iA- i(5 + CM - ^A^r - 



+-{B - C)Fo - djJ'Fo - 6kF2 + hjiJ'YFo + hjK^F^ + HkF^ 



(3) 



where J, Jj,, Jy, and are the molecular angular momentum operator and its projections on 
the molecular axes. The C, and A are the parameters of a rigid top which characterize the 
rotation around x, y, and z molecular axes, respectively (see. Fig. ||). The rest of parameters 
account for the centrifugal distortion effects [|l^]. In Eq. (|) the notation was used 



Fn — Jz {Jx 'Jv) ~^ {Jx '^y)Jz ■ 



(4) 



The Hamiltonian (^ can be diagonalized for each value of J and M in the basis of 
the symmetric-top quantum states \J,k,M > (—J < k < J), where J, k and M are the 
quantum numbers of the angular momentum, its projection on the molecular symmetry axis 
and on the laboratory quantization axis, respectively. The Hamiltonian (|^) has diagonal 
in k matrix elements due to the first two lines and matrix elements having \Ak\ = 2. The 
calculations of the rotational eigen states can be simplified if one uses the Wang basis 
with K = \k\: 



I"'"* =71 



\ao,p > 



1 

A 
1 + 



\a > +(-l)-'+^+P|a > 



0< K < J, 



\ao >; K = 0. 



(5) 



Here p = 0,1 and the sets of quantum numbers are a = {J, K, M}; \a >= {J, —K, M}; ao = 
{J,K = 0, M}. Depending on the parity of J, K and p, the states (^ generate 4 different 
irreducible representations of the molecular symmetry group D2 (Table 1). The molecular 



Hamiltonian is full symmetric (symmetry Ai). Consequently, the matrix elements between 
the states of different symmetry disappear. Thus diagonalization of the total Hamiltonian in 
the basis of (|) is reduced to the diagonalization of four independent submatrices, each for 
the states of particular symmetry. The rotational states of asymmetric top can be expanded 
over the basis states (|^) 

\P,P>=Y.^K\a,P>, (6) 

K 

where Ax stands for the expansion coefficients. The summation index, K, is given explicitly 
in (^, although depends on the other quantum numbers as well. All coefficients in 
the expansion (|^) are real numbers because the Hamiltonian is symmetric in the basis 
\a,p >. 

There are a few schemes for practical classification of the rotational states of asymmetric 
tops [|TB|. In this paper we will use the scheme which is somewhat better adapted to the 



description of nuclear spin isomers. We will designate the rotational states of asymmetric 
top by indicating p, J and prescribing the allowed K values to the eigenstates keeping both 
in ascending order. For example, the eigenstate having p = 0, J = 20, being the third in 
ascending order will be designated by /C = 4, because the allowed K in the expansion (||) 
are K = 0,2,4. . .20. (Note the difference between the two characters K and /C). It gives 
unambiguous notation of rotational states for each of the four species Ai, A2 (/C-even) and 
Bi, B2 (/C-odd). This classification establishes the connection with the prolate symmetric 
top for which IC = K. One should remember that physical meaning of the quantum number 
/C is clearly limited. To illustrate this we consider as an example the rotational states of Ai 
and A2 symmetry (/C-even) for J = 20. The upper panel in Fig. ^ shows that the energy of 
the rotational states is not determined solely by /C as it would be for a rigid symmetric top. 
The graph shows the difference in energy between rotational states Ai and A2 having the 
same /C-number but different quantum numbers p. The states Ai, /C = and A2, /C = are 
omitted from the graph because the latter state does not exist. As it is seen from the data 
the splitting is significant for low values of /C and rapidly disappears as /C increases. Note, 
that the Ai and A2 states for a rigid symmetric-top molecule would be degenerate for all JC. 

The low panel in Fig. ^ illustrates the property of the eigen state expansions over the basis 
states (^. In this panel the squared magnitude of the two Ak coefficients in each eigenstate 
of the symmetry Ai ( J = 20) is given. The first one ("/C-term" in the Fig. ^ is Aj^ having 
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K = JC. These coefficients appeared to be the biggest coefficients in each expansion. This 
adds some physical insight to the proposed classification scheme. The "second term" is the 
second biggest coefficient in the expansion. Again, at low /C values there are more than one 
significant terms in the expansion (|^). As /C grows, the contribution from just one term 
becomes predominant. Note that the same graph for a symmetric-top molecule would show 
only one significant term, Ak=ic, in the expansions of all eigenstates. 

IV. ORTHO AND PARA ISOMERS OF ETHYLENE 

The total molecular wave function is a product of spatial wave function and spin wave 
function. The nuclear spin states in ^"^00134 are of Ai and Bi species, having the statistical 
weights 20 and 12, respectively. The ortho spin states (symmetry Ai) have the two hydrogen 
pairs (H1-H2 and H3-H4) either both in triplet state, or both in singlet state. The para 
species (symmetry Bi) have one pair of protons in singlet state but the other pair in triplet 
state. Each symmetry operation of the group interchange in ^^CCH4 even number of 
protons. Consequently, the total wave function must be unchanged under all operations, 
thus belonging to the representation Ai. In order to have the total wave function being of 
species Ai one has to have the spatial wave function being of symmetry Ai and Bi as the 
spin wave functions are. This implies that the rotational states Ai and Bi should be only 
positive (even in parity) but the rotational states A2 and B2 should be only negative (odd 
in parity). 

Summarizing this discussion we can write the total states in the ^^CCH4 molecules. The 
ortho states can be presented as 

\fi>= \P,p> \ii2,i34.,ic >; ii2,k4. = t,s] /C - even. (7) 

Here ii2, ^34, and ic designate the spin states of the two pairs of protons and the nucleus 
^^C (Fig. 0). For ortho molecules the spin states of the proton pairs should be either both 
triplet, \t,t >, or both singlet states, \s,s >. 
The para states can be presented as 

1/^' >= 1/3', P > \i'i2, hi^^'c >; ^' - odd. (8) 

For para molecules one pair of proton is in singlet, but the other in triplet state, thus 
the proton spin states are \t',s' >, or \s',t' >. In (^, (§) and further we use unprimed 
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parameters for ortho species and primed parameters for para species; p = 0, 1 indicates 
positive, or negative sign of the state, respectively. Eqs. (|^, (|^) imply that the ortho-para 
state mixing in ethylene needs coupling of states having different parity of /C. 

The relative position of ortho and para states is important for the calculation of the 
conversion rate (see Eq. (]^)). In the case of a symmetric-top molecule the ortho-para energy 
gaps can be expressed as a polynomial of J and K. This gives "accidental" (isolated) 
resonances between the ortho and para states in symmetric top. In the case of an asymmetric 
top one has a phenomenon which can be called collapse of ortho and para states. It appears as 
a progressive decrease of the ortho-para energy gaps between the states of identical rotational 
momenta, J, as J increases. This is illustrated in Fig. ^, which shows the gaps between the 
^^CCH4 ortho states, /C = 0, and the para states, /C' = 1, both having positive sign {p = 0) 
and the same J. If such sequence of ortho-para level pairs would exist for symmetric tops it 
would dominate the conversion. In the case of asymmetric tops the situation is different. We 
will see below that this sequence of closed ortho-para level pairs do not produce a significant 
contribution to the isomer conversion in ^'^CCIl4. 

The Boltzmann factors Wo{<y) and Wp{a') in Eq. (|]) determines the relative population 
of rotational states in the ortho and para families. 

Pa = poWo{a); pa' = ppWp{a'), (9) 

where po and pp are the total densities of ortho and para molecules, respectively. The Eqs. (H) 
imply the equilibrium distributions inside the ortho and para subspaces. This is fulfilled with 
high accuracy even if the ratio po/ Pp is out of equilibrium because rotational relaxation is 
on many orders of magnitude faster than the ortho-para conversion. The partition functions 
for ortho and para molecules at room temperature (T=295 K) are found to be equal to 

Zortho = 2.66- 10^; Zpara = 1-60 ■ 10\ (10) 

In the calculation of these partition functions we took into account the degeneracy over M, 
nuclear spins, including also spin of nucleus ^^C, parity of states, as well as the restrictions 
imposed by quantum statistics. 
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V. ORTHO-PARA STATE MIXING 



In the present paper we will consider the ortho-para conversion induced by the two 
hyperfine interactions, viz., magnetic dipole-dipole interaction between the molecular nuclei 
(spin-spin interaction, Vss) and the nuclear spin-rotation interaction, V^^. Thus the total 
intramolecular perturbation able to mix the ortho and para states is 

V = Vss + VsR. (11) 
All matrix elements of V are diagonal in parity p. 



A. Nuclear spin-spin coupling 

The ortho-para conversion in molecules induced by nuclear spin-spin coupling were in- 
vestigated in (other references see in the review ^^). The spin-spin Hamiltonian for 
the two magnetic dipoles /xi and fi2 separated by the distance r has the form |TH 



^/f ^ = % - 3n.n, ; Pu = /ii/is/r^/^^^/^'^/i , (12) 

where i*-^-* and i*-^-* are the spin operators of the particles 1 and 2, respectively; n is the unit 
vector directed along r; i and j are the Cartesian indices. The second rank tensor T*^^^) in 
the Eq. (|12|) represents a spatial part of the spin-spin interaction. The second rank tensor 
acts on spin variables. 

The total spin-spin interaction in ^^CCH4 is composed from the interactions between all 
pairs of molecular nuclei. One can show that the spin-spin interactions between the protons 
1-3, 2-4, 1-2, and 3-4 have the spatial part which can mix the quantum states only if they 
have K numbers of the same parity. Consequently, these terms do not contribute to the 
ortho-para state mixing in ethylene. The complete nuclear spin-spin interaction in ^^CCH4 
able to mix the ortho and para states reads 

V,s = Vj?' + Vif + V<?' + V<?' + V<?' + Vfs'K (13) 

Here the upper indices refer to the hydrogen nuclei in the molecule 1. . . 4, and to the nucleus 
^^C, which has spin equal 1/2 (see Fig. Calculation of the spin-spin matrix elements can 
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be simplified by a few observations wliicli we discuss in more detail for the perturbations 
V^f^ and V^p. First, one can proof by applying proper symmetry operation that: 



<fi'\vg'^\fi>=<fi'\V^1,'>\f^> . 



(23), 



(14) 



Next, the two, out of four, matrix elements between the ortho and para spin states vanish. 



< t',s'\V^l^'>\s,s>=< s',t'\V^s'\s,s>= 0, 



(15) 



because in these two cases one has the matrix elements of a vector (spin operator) between 
the states both having zero spin. On the other hand, the remaining two matrix elements 
are equal to each other. 



<t\s'\vg'^\t,t>=< s',t'\vg'^\t,t>. 



(16) 



To summarize, one can conclude that the variety of matrix elements for the operators 
Vgg^^ and Vgf^ is reduced to one matrix element, e.g., < f3' ,p,t' , s'\Vgp\(3,p,t,t >. Obvi- 
ously, the carbon spin state, which is omitted in this expression for simplicity, should be 
unchanged in this matrix element. 

Further, we write this matrix element using an expansion over symmetric-top states (P) 



<P',p,t',s'\vg'^\P,p,t,t>=<t',s' 



A'j^,Ak < a',p\V^f^\a,p > 

K'K 



\tA > 



(17) 



This expression reduces the calculation of the spin-spin matrix elements in asymmetric tops 
to the calculation of symmetric-top matrix elements. Solution for the latter can be found in 
[P], which allows to express the strength of mixing in ethylene by Vss as 

/ 

Fss{a'\a) = {2f + l){2J+l)T^' 



K>0,q 




U±l)^ J' 2 J 

^ V2 -110 



Here q = ±1; (: : :) stands for the 3j-symbol and the notation was used 



(18) 



= mali'^l' + 2\Pci%T'\'' + 2\Pc3%T'\''; r = 46.5 kHz, 



-(CI) ,2 



(C3)|2. 



(19) 



where P-factors are equivalent to the similar factor in ([T^); T2^"'^ are the spherical compo- 
nents of the corresponding T-tensor calculated in the molecular frame. The numerical value 
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of T was found using the molecular structure from . The selection rules for the spin-spin 
mixing in ethylene read 



Ap = 0; |AJ|<2, (20) 
and in addition, parity of /C' and /C is opposite. 

B. Spin-rotation coupling 

Ortho-para conversion in molecules induced by spin-rotation coupling were studied in 



r9|, p0|] (more references can be found in [T^l)- Nuclear spin-rotation coupling in molecules 



is due to magnetic fields produced by the molecular electrical currents. The spin-rotation 
perturbation can be presented as ||2l|JT6|j20[| 

(21) 



VsR = ^ (E i^'^ • C^^^ • J + h.c^ ; « = 1, 2, 3, 4. 



For the spin-rotation perturbation relevant to the ortho-para mixing, the index i should 
refer only to the hydrogen nuclei. 

The calculation of the second rank spin-rotation tensor C is rather complicated problem 
which is not solved yet completely. We will use the following estimation of C. First we note 
that contribution to the spin-rotation coupling arising from the electric fields in the molecule 
has been shown to be very small compared to the part having "magnetic" origin and can 



be neglected [22|. Then, we split the C-tensor arising from magnetic fields of the moving 



charges into two parts, 

C = ^C + "C, (22) 

produced by the molecular electrons and nuclei, respectively. The nuclear contribution, "C, 
is obtained as a first order average in the vibrational ground state, and depends only on 
the nuclear coordinates. Contrary, the electron part, '^C, is a second order series expansion 
which involves the full electronic spectrum. Fortunately, xz and zx components which are 
responsible for the mixing of states having AK = 1 in C2H4 are vanishing by symmetry 



requirements 20 



Therefore, the C-tensor, effective in ortho-para conversion, can be written (in Hz) for 



the i-th as 22 
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bk = 2^pqk/chRl 



(23) 



where Rfc is the radial vector from the proton H*-*^ to the charge k] is the radial from the 
center of mass to the particle /c; qu are the nuclei' charges; B is the inverse matrix of inertia 
moment. B is a diagonal matrix having the elements B^^ = 58.6 GHz, Byy = 48.7 GHz, 
and B^z = 291.9 GHz. Index k runs here over all nuclei in the molecule except the proton i. 

The spherical components of the spin-rotation tensor of the rank I {I = 1, 2) for the i-th 
proton calculated in the molecular frame, Cj;^^, can be determined using Eq. (|23|). For the 



ethylene molecular structure from the Ref. [|18| and bare nuclei' charges these components 
are 



C^]} = 3.8 kHz; C{]1 = -2.9 kHz; 



Ci^l = -4.1 kHz; Ci;/ = 3.2 kHz. 



(3) 



(24) 



The difference between absolute values of C^^^ and C/^^ appears because of the shift of the 
molecular center of mass caused by the bigger mass of ^^C in comparison with ^^C. 

Similar to the previous section, one can reduce the calculation of the matrix elements 
of the spin-rotation coupling in asymmetric tops to the calculation of the symmetric-top 
matrix elements. The latter can be found in |T9| , pO| , p3| . Thus one has the the strength of 
mixing due to the spin-rotation coupling in ethylene. 



FsR{a'\a)=2{2J' + l){2J + l)Y^ 



1 + (-lY+p 



J2 A'j,^,AK^i;f,K + q\J,K) 

K>0,g 

2 



V2 



(25) 



Here g = ±l;i = l,3 denotes the hydrogen nuclei 1 and 3. Note, that the protons 2 and 4 
were also taken into account because they produce mixing equal to that of protons 1 and 3. 
In (^) the notation was used 

/ 



(0 

l,q 



V 



J' I J 

-K' q K 



X 



, . J' J l\ \ J J' I 

1 1 J 1 1 J' 



(26) 
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where {: : :} stands for the 6j-symbol; y{J) = ^Jj{J+ 1)(2J + 1). The selection rules for 
the ortho-para mixing by spin-rotation perturbation in ethylene read 



Ap = 0; |AJ| < 1. (27) 
And again, parity of /C' and /C is opposite. 

VI. CONVERSION RATES AND DISCUSSION 

The decoherence rate, F, is another important parameter for the isomer conversion by 
quantum relaxation (see Eq. (|])). This parameter is rather difficult to calculate or to 
measure. Due to its physical meaning the value of F is close to the rotational state population 
decay, which determines the line broadening in microwave rotational spectra. Unfortunately, 
the ethylene does not have pure rotational spectra because it has no permanent electric dipole 
moment. Or, to be precise, very small dipole moment in the case of ^^CCH4. 

In the present calculations we take the estimation of F from the line broadening in 
infrared spectra, for which the information is available. Recently the line broadening was 
accurately measured for the C2H4 molecules embedded in nitrogen One should not 



expect big difference for the line broadening in the isotope species C2H4 and CCH4. From 



the data ||2^ one can estimate the decoherence rate as 

T/P = 2-10^ s'VTorr. (28) 



This rate is 10 times smaller than the F determined for the isomer conversion in CH3F ||17|| . 
The decrease of F is because ethylene is not a polar molecule and thus has its rotational 
relaxation slower. Probably, the estimation (p8|) is too low for the case of pure ethylene gas 
because ethylene has polarizability larger than nitrogen. For the buffer gases having large 



polarizability, like Kr, or SFg, nearly three times bigger line broadening was found ||25|| . 
Nevertheless, we will use the estimation (^) which can be considered as a low limit for F. 

Now we are ready to calculate the conversion rate in ethylene. The final expression 
for the conversion rate, 7, is again the Eq. (§) having the strength of mixing, F{a'\a) = 
Fss{0''\0') + FsR{a'\a). The results of the calculations are given in the Table 2. From these 
data we note that there is no appreciable contribution to the conversion from the collapsing 
ortho-para level pairs presented in Fig. ^. Another observation is that the contribution from 
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the mixing of states having \A)C\ > 1 is also smalL This is the consequence of the fact that 
the ethylene molecule is rather close to a symmetric top. The main contributions to the 
conversion come from just two level pairs, both having the same quantum numbers (J',/C')- 
(J,/C), but different p. The most important ortho-para level pair (1,27,7)-(1,28,6) has the 
frequency gap ^1.0 GHz. This level pair is indicated in Fig. ^ from which one can see that 
this pair of states is situated at rather high energies. Nearly 10% of the conversion rate is 
due to the pair (0,27,7)-(0,28,6). This level pair has the same strength of mixing as the first 
pair but three times bigger level splitting, u. The properties of the states (1, J = 28, /C = 6) 
and (1, J = 27, /C = 7) are illustrated in Fig. ^ which presents the expansion coefficients. Ax, 
for these states. One can note that the biggest coefficient in both expansions has the value 
of K equal to the value of /C chosen for the designation of the eigenstate in our classification 
scheme. 

The total conversion rate in ^'^CCH4 which is the sum over all ortho-para level pairs 
having J' and J up to 40, was found equal 

7/P = 2.7 ■ 10"^ s^VTorr. (29) 

This rate is close to the experimental value (5.2 ± 0.8) ■ 10^'^ s^^/Torr 0. 

The calculations of the conversion rates were repeated using less accurate molecular 
parameters from [|T^]. We obtained essentially the same results. The conversion is again 
determined by the same two level pairs and the value of the total conversion rate is 3.5 ■ 
10~^ s~^/Torr, which is close to the more precise value (pQ]). 

Our model allows to determine the pressure dependence of the conversion rate in ^^CCH4 
at the conditions of the experiment [^. By comparing the level splitting of the most impor- 
tant ortho-para level pair (1,27,7)-(1,28,6), which is ~1.0 GHz, with the decoherence rate 
( PSD one can conclude that the case corresponds to the limit F ^ cj. From Eq. (^ one 
can see that in this limit linear grows of the conversion rate, 7, versus gas pressure should 
take place. (Note that F is proportional to the gas pressure.) This dependence was indeed 
observed in the experiment 0]. 

Magnitude of the calculated conversion rate relies heavily on the determination of the 
ortho-para level splitting. The gap of the most important level pair appeared to be rather 
large, ~1.0 GHz, and is unlikely that it has been determined with significant error. This can 
be proven by comparing this splitting with the one calculated using less accurate molecular 
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parameters [0. These parameters give the gap 0.914 GHz, which is close to the value given 
in the Table 2. 

The two sources of the ortho-para mixing have been analysed. The spin-spin interaction 
between the molecular nuclei was possible to calculate rather accurately. The uncertainty of 
this perturbation is mainly due to small errors in the knowledge of the molecular structure. 
On the other hand, our analysis has shown that the spin-spin coupling contributes less 
than 1% to the ethylene conversion and thus negligible. There are a few reasons for this 
reduction. First, the spin-spin tensor has relatively small value due to the big distances 
between interacting nuclei and small magnitude of the magnetic moment of the carbon 
nucleus. Second, due to high symmetry of the molecular structure interactions between 
many pairs of protons do not contribute to the ortho-para mixing. Moreover, there are 
no close ortho-para level pairs at small J, where spin-spin coupling can compete with the 
spin-rotation coupling. 

We have taken into account also the spin-rotation perturbation. The magnitude of the 
spin-rotation tensor was estimated and it was demonstrated that the main contribution 
to the ethylene conversion comes from the spin-rotation mixing. Strong contribution from 
the spin-rotation coupling is due to the fact that the strength of the mixing Fsr grows 
as in comparison with the slower (J^) grows of the spin-spin strength of mixing Fss- 
This makes the "accidental" resonance at big J (1,27,7)-(1,28,6) being so important. An 
extra cause of the efficiency of the spin-rotation coupling is that the spin-rotation tensors 
of first and second rank contribute both to the conversion. On the other hand, not all close 
ortho-para level pairs contribute to the conversion. In the particular case of the collapsing 
ortho-para levels (Fig. |^) there is no strong mixing mainly because of interferences between 
many symmetric-top components with different quantum numbers K, and also because for 
these J' = J pairs the first rank spin-rotation tensor does not contribute. 

Error in the estimation of F introduces systematic uncertainty into the final result. But 
it does not undermine the model as a whole. First, one can turn around the approach and 
determine which decoherence rate F should be chosen in order to fit the experimental data. In 
order to reproduce the experimental rate of conversion one should take the decoherence rate, 
r/P = 3.9 • 10^ s~^/Torr, which is less than two times different from our rough estimation 
Of course, such approach will be justified only when the spin-rotation perturbation 
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in ethylene will be known from an independent source. There is another point which is 
worth to mention. Recently it was proposed a new approach to the isomer conversion in 
which the knowledge of F is not necessary at all |2^. It is based on a fast linear sweep 



of the molecular levels through the ortho-para resonance. Nuclear spin conversion in such 
experimental arrangement does not depend on F and is dependent on the strength of ortho- 
para mixing and the population of the mixed states. 

Finally, it seems interesting to emphasize the difference in the conversion rates produced 
by the two pairs of identical rotational quantum numbers (27,7)-(28,6) but different in 
parity. This peculiarity arises from the asymmetric top properties, where energy of states 
depend on parity. This in turn has dynamical consequence that the spin isomers equilibrate 
first in one parity manifold (here the odd one). 

VII. CONCLUSIONS 

The theoretical model for the nuclear spin conversion in ethylene (^^CCH4) has been 
developed. For the first time a theory of the spin isomer conversion in asymmetric tops was 
possible to compare with experimental data on the isomer conversion. We have found that 
the two experimental results the magnitude of the ethylene (^^CCH4) isomer conversion 
rate and its pressure dependence, are consistent with the spin conversion governed by quan- 
tum relaxation. The ortho-para state mixing is performed in this molecule mainly by the 
coupling between the protons' spins and the molecular rotation. We have identified also the 
two pairs of ortho-para states which are almost completely determine the spin conversion in 
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Table 1. The character table for the symmetry groups C2v, C2,,(M), D2 and classification 
of the basis states (|^). 



C2.(M) 


E (12)(34) E* (12)(34)* 


ortho 


ortho para 




E 


C2 








D2 


E 


K 






K-eyen 


K=0 K-odd 


Ai 


1 


1 


1 


1 


p=0 


J, p-even, - 


B2 


1 


-1 


-1 


1 




- p=l 


A2 


1 


1 


-1 


-1 


p=l 


J, p-odd, - 


Bi 


1 


-1 


1 


-1 




p=0 
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Table 2. The most important ortho-para levels and their contributions to the spin con- 
version in ethylene. In calculations the molecular parameters from WM were used. The total 
rate combines the contribution from all ortho-para level pairs having J < 40. 



Level pair 


Energy 


LO /2tx Fss 


F'sR 


l/P 


p', J', IC'-p, J, K. 


(cm^-*^) 


(MHz) (MHz2) 


(MHz2) 


(10-^ s-VTorr) 


1,27,7-1,28,6 


871.53 


1006.7 4.2-10-2 


4.4 


2.49 


0,27,7-0,28,6 


871.69 


-3435.5 4.2-10-2 


4.4 


0.21 


0,4,3-0,6,2 


53.72 


-3180.6 1.1-10-3 





3.4-10-3 


1,26,3-1,24,6 


680.68 


944.3 3.4-10-^ 





5.6-10-^ 


0,20,1-0,20,0 


356.02 


1894.4 7.3-10-^ 


2.4-10-6 


4.8-10-6 


Total rate 








2.70 
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ortho para 

FIG. 1. The ortho and para states of the ^^CCJ?4 molecules. The bent lines indicate the 
rotational relaxation inside the two subspaces. Va'a refers to the intramolecular mixing of the ortho 
and para states. The indicated pair of states is the most important one for the spin conversion in 
^^CCH4. 
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FIG. 2. Numbering of atoms in ^^CCH^ and orientation of the molecular system of coordinates. 
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FIG. 3. Properties of asymmetric-top rotational states. Upper panel gives splitting between 
the states (0, J = 20, IC) and (1, J = 20, /Cj. Low panel gives squared values of the two terms in 
the expansion for the state (0,J = 20, IC). "IC-term" is A\^^. "Second term" is the second 
biggest coefficient in each expansion. 
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FIG. 4. Collapse of the ortho and para states in ethylene. The figure shows the energy gaps 
between the ortho states (0, J, IC = 0) and the para states (0, J, /C = 1 j. 



24 



0.8 



g 0.4 



0.0 



-0.4 



o 



17 



o 



b#o#o#o#0#o#o#o#c 







10 



20 



K 



FIG. 5. The expansion coefficients, Ak, for the states most important for the spin conversion 
in ^^CCH4. (o)-ortho state (1, J = 28,/C = 6); (•)-para state (I, J = 27,IC = 7). 
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